Constraints on Minimal SUSY models with warm dark matter neutralinos 
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If the energy density of the Universe before nucleosynthesis is dominated by a scalar field cf> that 
decays and reheats the plasma to a low reheating temperature Trh , neutralinos may be warm dark 
matter particles. We study this possibility and derive the conditions on the production mechanism 
and on the supersymmetric spectrum for which it is viable. Large values of the parameter and 
of the slepton masses are characteristic features of these models. We compute the expected direct 
detection cross sections and point out that Split-SUSY provides a natural framework for neutralino 
warm dark matter. 
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Thermally produced neutralinos are typical and well- 
motivated cold dark matter candidates (see e.g. pj). 
They are produced by scatterings in the thermal bath, 
then reach equilibrium and finally decouple when non- 
relativistic. Non-thermally produced neutralinos, on the 
contrary, are usually produced as relativistic final states 
in the decay of heavy particles and may never reach chem- 
ical or kinetic equilibrium. As a result, they are not 
necessarily cold. Indeed, if they manage to keep most of 
their initial energy they might behave as warm dark mat- 
ter suppressing the evolution of small-scale structures in 
the Universe. 

Although consistent with the observations of the large 
scale structure of the Universe and the cosmic microwave 
background radiation anisotropies, cold dark matter 
models seem to have problems at galactic scales. They 
not only tend to form cuspy structures in the halo density 
profile [2| but also predict a large overabundance of small 
halos near galaxies such as our own Q . Warm dark mat- 
ter, with its larger free-streaming scale, may solve these 
problems while maintaining the celebrated success of cold 
dark matter models at large scales 

Neutralinos as warm dark matter candidates were ini- 
tially discussed in Refs. 0.0 In this paper we will ex- 
tend, in several ways, the analysis presented in those ref- 
erences. Ref. concentrated on the cosmological aspects 
of this possibility, while we study a particular particle 
model. Reheating temperatures smaller than 5 MeV or 
2 MeV as those found in Ref. are hardly compatible 
with the standard cosmological scenario, which require 
T RH > 4 MeV 0. We, instead take T RH ~ 10 MeV as a 
characteristic value and find the conditions on the initial 
energy and on the supersymmetric spectrum for which 
the neutralino is a warm dark matter particle. In doing 
so, we ignore, following Ref. [8j, the naturalness argu- 
ment giving up the idea that SUSY stabilizes the weak 
scale, notion that is the basis for Split-SUSY models. We 
also discuss the non-thermal production of neutralinos, as 
well as the implications for direct dark matter searches 
of neutralino warm dark matter. 

We concentrate on non-standard cosmological models 
(see for example Ref. |9() in which the late decay of 
a scalar field <j> reheats the Universe to a low reheat- 
ing temperature Tru, smaller than the standard neu- 



tralino freeze out temperature. Such scalar fields are 
common in superstring models where they appear as 
moduli fields. These fields get mass at the low energy 
supersymmetry breaking scale, typically of the order of 
10 2 - 10 3 TeV. The decay of (j> into radiation increases the 
entropy, diluting the neutralino number density. The de- 
cay of <p into supersymmetric particles, which eventually 
decay into neutralinos, increases the neutralino number 
density. We denote by b the net number of neutrali- 
nos produced on average per <f> decay. The number b 
is highly model-dependent, so is the <fi field mass m^. 
They are determined by the physics of the hidden sec- 
tor, by the mechanism of supersymmetry breaking, and 
in superstring-inspired models by the compactification 
mechanism [a llfj) . Thus, here we consider b and as 
free parameters. 

To account for the dark matter of the Universe, the 
neutralino relic density must be in agreement with the 
observed dark matter density. In low Trh cosmologi- 
cal models essentially all neutralinos can have the dark 
matter density provided the right combination of the fol- 
lowing two parameters can be achieved in the high energy 
theory: the reheating temperature Trh, and the ratio of 
the number of neutralinos produced per <f> decay over the 
(j) field mass, i.e. b/m^ 0, EJ. We will find later the 
values of Trh and b/m^ for which the neutralinos we 
are interested in have the right dark matter density. 

A crucial quantity that distinguishes warm from cold 
dark matter is the free-streaming length at matter- 
radiation equality \fs, which depends on the param- 
eter r x — a(t)p x (t)/m x (see for example Refs. 0, Q 
and references therein). This parameter would be the 
present characterictic speed of neutralinos of mass m x , if 
their momentum p x only redshifted from neutralino pro- 
duction onwards (a(t) is the scale factor with ao = 1). 
During the cosmic evolution r x is constant. Structures 
smaller than Xps are damped because neutralinos can 
freely flow out of them. N-body simulations have shown 
that to explain the lack of substructure in the local group, 
Xfs should be of order 0.1 Mpc 0, thus r x ~ 10" 7 . 

The parameter r x also determines (for neutralinos 
which are relativistic at production 12]) the neutralino 
phase-space density Q. This is defined as Q = p/{v 2 ) 3 / 2 
where p is the neutralino energy density and (v 2 ) is the 
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mean square value of the particle velocity. In the ab- 
sence of dissipation, the coarse-grained phase-space den- 
sity (the quantity that can actually be observed) can 
only decrease from its primordial value. The observa- 
tion of dwarf-spheroidal galaxies places a lower bound on 
Q which translates into an upper bound on r x of about 
2.5xl0~ 7 In what follows we impose the condition 



r x — 10 7 to our model. 



Neutralinos are produced in </> decays, thus their aver- 
age initial energy is Ej = m^/N, where N is a, number 
which depends on the production spectrum. We expect 
N to be of order one and require Ej » so that neu- 
tralinos are relativistic at production. Thus, p x ~ Ei at 
the moment of <j) decay. Assuming an instantaneous <f> 
decay at Trh, with no subsequent entropy production, 
the scale factor at decay is a = Tq/Trh, where To is the 
present photon temperature, and the parameter r x in our 
model is r x — (T Ej) / (T RH m x ) i.e. 
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Thus, the condition r x = 10 7 fixes in terms of the 
reheating temperature and the neutralino mass 



m = lO^TeV 
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In the estimation of r x we have assumed that the neu- 
tralinos do not lose their energy in scattering processes 
with the thermal bath. To ensure this condition we will 
simply require that the interactions of neutralinos with 
the particles present in the plasma, e ,v and 7, are out of 
equilibrium. Neutralino interactions are determined by 
the neutralino composition in terms of gauge eigenstates. 
The ligthest neutralino (x) is a linear superposition of 
bino (B°), wino (W°) and higgsino (H°) states, 
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NuB° + N 12 W° + N 13 H° + N 1A M . 
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According to the dominant term in Eq. © , x is classified 
as bino-like, wino-like, or higgsino-like. 

Wino and higgsino-like neutralinos are always accom- 
panied by a chargino state (x ± ) with a mass difference 
Am = m x ± — m x ~ sin 2 (8 'w)(Mz/ni x ) 2 m x for higgsi- 
nos and much smaller for winos. When m x Am < EjT 
relativistic neutralinos scatter inelastically (e.g. %e~ — > 
X~v) besides scattering elastically (e.g. x v ~~ > X v ) on 
electrons and neutrinos. Since it is not possible to simul- 
taneously suppress the elastic and the inelastic scattering 
rates, wino and higgsino-like neutralinos easily lose their 
initial energy. In Ref. Q was indeed pointed-out that 
a wino-like neutralino cannot be warm dark matter and 
that a higgsino-like neutralino would only be viable if 
Trh < 2 MeV were allowed. Bino-like neutralinos, on 
the other hand, couple to electrons and neutrinos mainly 
through slepton exchange and typically have small inter- 
action cross sections. In general, there are no charginos 



close in mass to bino-like neutralinos and inelastic scat- 
terings can be suppressed. We will therefore focus on 
bino-like neutralinos in the following. 

In the approximation > 6TEj, which is guaranteed 
by the lower bound found below, the cross section for 
elastic bino-lepton (x-£) scattering mediated by a slepton 
of mass mz is estimated to be 



g 4 EtjEjf 
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Here g is the weak hypercharge coupling constant, En ~ 
3T is the energy of the incoming thermal-bath neutrino 
or electron, and Eg / is the energy of the outgoing lepton. 
We require that interactions are out-of-equilibrium, 
i.e. that the interaction rate is smaller than the expan- 
sion rate of the Universe, T ~ nav < H, where n is the 
thermal neutrino or electron number density and H is 
the Hubble parameter. Since the interaction rate goes as 
T 5 and H as T 2 , if the interactions are out-of-equilibrium 
at Trh they will also be out of equilibrium at any later 
temperature T < Trh- Replacing Egf in Eq. 0] by its 
maximum value Ej = m^/N, we obtain that neutralino 
interactions are out of equilibrium at Trh (and subse- 
quently) if 



RH 



0.1 ( m * V 

ViV10 3 TeVy V 10 MeV 



T 



2 'AO TeV\ 4 



< 1. 



(•5) 



Sleptons therefore must be heavy. 

Due to the structure of the neutralino mass matrix, 
bino-like neutralinos always have a small higgsino com- 
ponent. As a result, they also couple to fermions through 
Z exchange. If sleptons are heavy, the Z-exchange dia- 
gram could give the dominant contribution to the scat- 
tering rate. In this case the neutralino cross section is 



ov < tan 4 
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where 9yv is the weak mixing angle and tan/3 = v-z/v\ is 
the ratio between the vacuum expectation values of the 
two Higgs doublets. Eq JBJ holds for m| < 6TE1 and 
m x ^ 6TEj. The latter condition, after using Eq. [21 
becomes T RH /10MeV > 1.3(m x /100GeV). The condi- 
tion r < H implies for the interaction in Eq. H3 (again 
replacing Ei>f by its maximum possible value, i.e. Ex) 
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If m 2 > 6T£ 7 (i.e. T RH /10MeV < 1.3(m x /100GeV), 
the right-hand side of Eq. El and the left-hand side of 
Eq. [7| must be multiplied by the factor m x /3TRHEi, 
which yields 
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In any event, a large \x parameter value is required. 

Replacing Eq. |21 into Eqs. and [7| and directly from 
Eq. 03 we obtain bounds of the slepton masses and the 
parameter fi which depend only on the neutralino mass 
and the reheating temperature, 

•"'•>—(To^v) 1,4 (sfv) 3 ' 4 - < 9 > 
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from Eq. 0or 
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when Eq. [HI holds instead. 

Using the cosmological lower bound Trh > 4 MeV we 
obtain m E > 9 TeV (mx/lOOGeV) 1 / 4 and fx > 2.0 TeV 
(rnx/lOOGeV) 1 / 4 or p, > 2 TeV (mx/lOOGeV) 1 / 2 . 

Let us now return to the density requirement. In gen- 
eral there are four different ways in which the neutralino 
relic density flh 2 depends on Trh, in the low reheat- 
ing cosmologies we are considering. The four cases have 
thermal or non-thermal dominant neutralino production, 
with or without chemical equilibrium [ToLIIlj . The case of 
interest here is that of a very low reheating temperature 
and high standard neutralino relic density (corresponding 
to small neutralino cross sections) for which the produc- 
tion is non-thermal without chemical equilibrium [Tcj . In 
this case neutralinos are produced in the decay of the <j> 
field into supersymmetric particles, with subsequent fast 
decay of all other supersymmetric particles into neutrali- 
nos, and the production is not compensated by annihila- 
tion. This implies 

thJfL ~ o 10 3 6 f 103TeV } ( m x ) ( T ™ \ 
0.11 V m 4> / V10 2 GeV/ V 10MeV / 

This equation fixes the value of b/m^ required for neu- 
tralinos to have the dark matter density, fld m h 2 = 0.11 1 , 

m ^ TeV J VlOOGeV/ \ T RH J ' [ ' 

Using the program described in Ref. we have fol- 
lowed the evolution of the relic density of some neutrali- 
nos warm dark matter as function of Trh- This is shown 
in Fig. 2] We use a Minimal Supersymmetric Standard 



1 Cl dln h 2 = 0.109io'So6 was obtained for a ACDM model with 
scale-invariant primordial perturbation spectrum through a 
global fit of cosmic microwave background, supernovae, and large 
scale structure data Il4l 
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FIG. 1: The neutralino relic density as a function of Trh in 
four different models compatible with warm dark matter. In 
all of them M 2 ,M 3 ,mA = 1 TeV, tan/3 = 10, A t ,A b = 0. 
In A, m x ~ Mi = 100 GeV, /x = 10 TeV, m f ~ = 50 TeV, 
and b/m^ = 4 x 10" 7 TeV" 1 . In B, m x ~ Mi = 100 GeV, 
H,m~ f = 100 TeV, and bjra^ = 4 x 10" 7 TeV" 1 . In C, m x ~ 
Mi = 300 GeV, fi = 10 TeV, mj = 50 TeV, and b/m^ = 1.3 x 
10" 7 TeV" 1 . In D, m x ~ Mi = 300 GeV, (i,mj = 100 TeV, 
and b/m^ = 1.3 x 10" 7 TeV" 1 . 



Model defined in terms of the parameter set M3, Mi, 
Mi, rriA, fx, tan/3, mo, At, and A\,. Here Mj are the 
three gaugino masses, tua is the mass of the pseudoscalar 
higgs boson, and tan j3 denotes the ratio V2/V1. The soft 
breaking scalar masses are defined through the simplify- 
ing ansatz Mq = M\j = Mr> = Me = Ml = niQ whereas 
the trilinear couplings are given by Ajj = diag(0, 0, A t ), 
Ad = diag(0, 0, Ab), and Ae = 0. All these parameters 
are defined at the weak scale. The models presented in 
Fig. [U have m x ~ Mi = 100 GeV (solid lines) or 300 
GeV (dashed lines) and other parameters as mentioned 
in the figure caption. 

Fig. ^ confirms that at low reheating temperatures 
the neutralino production is purely non-thermal, inde- 
pendently of the supersymmetric spectrum, and follows 
Eq. (|12f> . Thermal effects are relevant only for Trh > 
1 GeV. The neutralino relic density shown in the figure 
correspond to bjm^ = 4 x 10" 7 TeV" 1 for m x — M\ 

1 ill) GeV (solid lines), and to b/m^ = 1.3 x 10" 7 TeV" 1 
for m x = Mi = 300 GeV. 

Combining Eas.mand ll2l the requirement of neutrali- 
nos to account for the whole of the dark matter and to be 
warm dark matter, fixes b and (see Eq.[2J) separately, 

^ 5 io=(A)(^). ,») 

That is, fixing fl x h 2 — 0.11 and r x — 10" 7 requires b = 

2 10" 4 N. 

Neutralinos, therefore, can be warm dark matter par- 
ticles only if slepton masses and \x are both large, above 
about 10 and 2 TeV respectively. Supersymmetric mod- 
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FIG. 2: Spin-independent neutralino-proton cross section for 
warm dark mater neutralinos as function of the squark mass, 
for different values of niA, and for Mi = m x = lOOGeV, 
Ma = M 3 =lTeV, fi = lOTeV, slepton masses of 50TeV and 
tan/3 = 10. The horizontal dotted line shows the reach of 
Zeplin IV-MAX and Xenon IT. 

els with heavy sfermions do exist in the literature. In 
Split-SUSY models 0, E3> E| , all scalar superpartners 
are usually heavy, with a mass scale that, in principle, 
could go up to the GUT or the Planck scale. These mod- 
els solve a number of phenomenological problems associ- 
ated with ordinary supersymmetry: they alleviate proton 
decay, increase the mass of the Higgs boson, and mini- 
mize the flavor and CP problems. Thus, a Split-SUSY 
model with a bino-like LSP and fi > 10 TeV would be 
a perfectly viable framework for neutralino warm dark 
matter. 

In standard cosmological models, a bino-like neutralino 
is allowed in Split-SUSY models only if it is almost degen- 
erate in mass with other neutralino and chargino states 
(Mi ~ M2 or Mi ~ u ) an d coannihilation effects deter- 
mine its relic density |16| . These constraints, however, 
do not hold in our model. By assuming a non-standard 
cosmology which allows for non-thermal production of 
neutralinos we have effectively enlarged the viable pa- 
rameter space of Split-SUSY models. 

In Split-SUSY models all scalar superpartners are usu- 
ally heavy but our conditions require only that the slep- 
tons be heavy, not necessarily the squarks. A model with 
light squarks would certainly be appealing from a phe- 
nomenological point of view. The most general low en- 



ergy spectrum of masses below a TeV of a model compat- 
ible with neutralino warm dark matter would consist of: 
the lightest neutralino, the next-to-lightest neutralino, 
the ligthest chargino, the three Higgs bosons, the gluino 
and the squarks. Such low energy spectrum is sufficiently 
rich to give significant signals at the LHC, particularly 
for light gluinos and squarks. Besides, the absence of 
sleptons in the low energy spectrum would be a crucial 
test for models with neutralino warm dark matter. 

The prospects for direct detection of bino-like neu- 
tralinos depend on the neutralino-nucleon cross section, 
which in turn is determined by the diagrams mediated by 
the heavy neutral Higgs boson (H) and squarks. We show 
in Fig.Elthe spin-independent neutralino-proton scatter- 
ing cross section as a function of the squark mass for 
three different models compatible with neutralino warm 
dark matter, each with a different value of tua, the pseu- 
doscalar Higgs boson mass. The mass of heavy boson H 
is close to rriA- The lines are not continued to smaller 
squark masses because they enter into a region incom- 
patible with the higgs mass bound 0] or b — > S7 [Tfll |. 
The figure shows that the expected neutralino-nucleon 
scattering cross sections in neutralino warm dark mat- 
ter models are small, between 10 -46 cm? and 10 -49 cm? . 
Hence, the expected spin-independent neutralino-proton 
cross sections are below the reach of the largest dark mat- 
ter detectors envisioned at present, i.e. Zeplin IV-MAX 
and Xenon 1 Ton poj . 

In this letter we studied non-thermally produced neu- 
tralinos that are viable warm dark matter candidates, 
in low reheating temperature cosmological scenarios. 
We have shown that models compatible with neutralino 
warm dark matter satisfy a number of non-trivial require- 
ments, some of them testable: the reheating temperature 
is low, the mass of the decaying scalar field which re- 
heats the Universe and the number of neutralinos it pro- 
duces per decay are constrained, the neutralino is bino- 
like, slepton masses and the mu parameter are both large 
and the neutralino interacts too weakly to be observed 
in direct dark matter experiments in the near future. To 
validate the idea of neutralino warm dark matter, there- 
fore, cosmological observations as well as accelerator and 
dark matter searches will be essential. 
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